Il Paesaggio archeologico dell’Asia Occidentale Antica

Studio e analisi dei sistemi di indagine remote sensing, morfologia e topografia

Andrea Titolo

Università di Torino

Muoversi nelle slide

  • F tutto schermo
  • → , ↓ , SPACE , N, L , J Slide successiva
  • ← , ↑ , P , H , K Slide precedente
  • Shift + ←/↑/→/↓ Vai alla prima/ultima slide
  • ESC, O Panoramica delle slide
  • M Menu

Introduzione ai GIS - Georeferire una mappa

Introduzione ai GIS - Georeferire una mappa


Inizieremo ad affrontare la questione della georeferenziazione di mappe ed immagini satellitari, lavorando su diverse tipologie di mappe e situazioni.

  • Principi di georeferenziazione e digitalizzazione di una mappa
  • Georeferire una mappa topografica con coordinate precise
  • Georeferire una mappa topografica non precisa
  • Georeferire una mappa geomorfologica usando punti già esistenti
  • Georeferire una foto satellitare


Nel frattempo, aprire QGIS!

Link alle slides:

Una volta caricate dal link, le slide dovrebbero seguire la presentazione in automatico

Georeferire una mappa – la teoria

Cosa vuol dire georeferire?

  • Assegnare coordinate geografiche ad un’immagine che non ne possiede, tramite il processo di gereferenziazione (o georettificazione, diverso da ortorettificazione)
  • Molte mappe (storiche, ma anche ad es. geologiche) non possiedono informazioni locazionali.
  • Quando vengono caricate in GIS, il software considera la grandezza dei pixel come 1m quadrato e e posizionerà la mappa scansionata alle coordinate 0,0 nell’oceano atlantico (NULL ISLAND)
  • Per evitarlo bisogna assegnare delle coordinate geografiche ai pixel
  • Il fine ultimo è quello di digitalizzare informazioni (cartacee) all’interno di una piattaforma GIS, e farle dialogare con altre informazioni spaziali.

Quali input vengono utilizzati?

Se si scansiona una mappa personalmente, assicurarsi di scansionare le immagini a 300dpi per una risoluzione ottimale. Altrimenti, ci si deve accontentare di quello che si ha.

Le mappe possono essere svariate:

  • Mappe storiche (a diverse scale)
  • Mappe di distribuzione
  • Mappe geologiche
  • Mappe geomorfologiche
  • Mappe topografiche
  • Foto aeree/drone
  • Foto satellitari

Come si fa nella pratica?

  • Si sostituiscono le coordinate dei pixel con valori geografici cliccando su un luogo visibile sia nella mappa da georeferire, che nella mappa di base in QGIS.
  • Questi punti di riferimento che si prendono sono definiti Ground Control Points (GCPs)
  • A seconda della quantità di informazioni in una mappa, la gereferenziazione può essere più o meno semplice
  • Una mappa con elementi chiaramente definiti come edifici, strade, torri, ecc. è molto più facile da georeferire rispetto ad una mappa geologica

Quale mappa di base si usa?

  • Ovviamente, la mappa che si utilizza come riferimento deve essere già georeferita
  • È possibile utilizzare altre mappe storiche, topografiche, o simili, di cui si è certi dell’accuratezza della georeferenziazione
  • Il più delle volte si utilizzano immagini satellitari a media risoluzione (ad es. Landsat, Sentinel) per aree regionali. o immagini ad alta risoluzione per aree più piccole (SPOT, Worldview)
  • Altrimenti, si utilizza la mappa di base disponibile in QGIS (Google, Bing, Esri, OpenStreetMap)
  • L’importante è utilizzare la stessa mappa, e non prendere punti su una o sull’altra (a causa di piccoli problemi di distorsione delle immagini)
  • Se si hanno più immagini dello stesso tipo da georeferire, utilizzare sempre la stessa fonte come mappa di base per la georeferenziazione

Quali elementi sono utili come Ground Control Points?

  • I GCPs possono anche essere raccolti sul campo
  • La scelta dei GCPs è fondamentale, perchè permettono alla mappa di posizionarsi, migliori GCPs permettono una migliore accuratezza
  • Per georeferire una mappa bastano anche solo 4 GCPs, tuttavia alcuni casi estremi ne richiedono tra 200 e 300.
  • Principi per buoni GCPs:
  • Ben visibili su entrambe le immagini
  • Buon contrasto o contorni ben definiti rispetto agli elementi circostanti
  • Oggetti che non si muovono nel tempo

Quali elementi sono utili come Ground Control Points?

  • Esempi di buoni GCPs
    • Edifici.
    • Evitare tetti degli edifici (se possibile) e preferire gli angoli
    • Angoli o incroci di strade
    • Angoli di altri elementi
    • Punti di riferimento: Ponti, Torri, Monumenti
    • Siti archeologici
    • Con cautela: tratti di fiume, elementi del paesaggio, ecc. (controllare che siano presenti su entrambe le immagini)

Cosa vuol dire digitalizzare?

  • La digitalizzazione è il processo tramite cui gli elementi di una mappa cartacea sono trasferiti in formato digitale (shapefile, geopackage, ecc.)
  • La digitalizzazione è il passaggio successivo alla georeferenziazione (non sempre).

Principi di digitalizzazione

  • La digitalizzazione (come la gereferenziazione) comporta sempre delle imprecisioni, e sarà impossibile ottenere un dataset corretto al 100%
  • Tuttavia, si può lavorare per ridurre al minimo le imprecisioni:
  • Usare lo stesso sistema di riferimento per mappe di base e elementi digitalizzati
  • Considerare la scala e la risoluzione della mappa di base
  • Scegliere la corretta geometria per i singoli elementi
  • Considerare i dati già a nostra disposizione
  • Considerare la qualità dei dati raccolti e la compatibilità tra i diversi dati
  • Provengono da raccolte sul campo? sono stati presi da qualche fonte online? in che modo sono stati raccolti originariamente

Georeferire una mappa in QGIS

Georeferire una mappa in QGIS

Prima di tutto, per facilitare l’orientamento, carichiamo in QGIS un layer che ci indichi dove si andranno a posizionare le nostre mappe (questo layer ci aiuta in questa situazione, ma in un progetto personale si dovrebbe ovviamente già essere a conoscenza dell’area coperta dalla mappa che si è scaricata).

Apriamo quindi il nostro gestore delle sorgenti dati e andiamo su Vettore. Sotto sorgente apriamo la ricerca e scegliamo lo shapefile nominato bbox_mappe.shp.

Georeferire una mappa in QGIS

Una volta scelto, clicchiamo su Aggiungi.

Georeferire una mappa topografica con coordinate precise

Georeferire una mappa topografica con coordinate precise

Una volta caricato il file, apriamo la tabella attributi del layer. Scegliamo (cliccando con il numero corrispondente alla riga della tabella) la geometria con nome Bar-Zertal2021 e clicchiamo su Zoom alla selezione (il pulsante con la lente d’ingrandimento).

Questo sposterà la visuale della nostra mappa sulla zona coperta dalla geometria selezionata.

Georeferire una mappa topografica con coordinate precise

Una volta fatto, ci troveremo nel Levante Meridionale. Procederemo a georeferire una mappa di distribuzione di siti archeologici con base una mappa topografica. La mappa presenta già una griglia di coordinate che in questo caso ci aiuteranno a georeferire velocemente l’immagine.

Georeferire una mappa topografica con coordinate precise

Procediamo quindi ad aprire lo strumento per la georeferenziazione. Avremo bisogno dello strumento “Georeferenziatore Raster”. Il Georeferenziatore è uno strumento nativo di QGIS (viene installato assieme al programma). Originariamente un plugin, nelle versioni meno recenti di QGIS deve essere attivato (se non lo è già) sotto Plugin → Gestisci ed Installa Plugin → Installati e attivarlo.

Per accedere al georeferenziatore, clicchiamo su Layer → Georeferenziatore (oppure sotto Raster → Georeferenziatore)

Georeferire una mappa topografica con coordinate precise

Ci si aprirà una nuova finestra, il Georeferenziatore. Qui potremo caricare la nostra immagine, i nostri GCPs e modificare le impostazioni di georeferenziazione (che vedremo in seguito). Per prima cosa, carichiamo l’immagine da georeferire, clicchiamo sul pulsante Apri Raster in alto a sinistra

Georeferire una mappa topografica con coordinate precise

Selezioniamo il file denominato “Bar-Zertal2021.png” e clicchiamo su apri. Questa mappa proviene da una di molte pubblicazioni della serie “Manasseh Hill Survey” condotte da Adam Zertal a partire dagli anni ’90.

Georeferire una mappa topografica con coordinate precise

Una volta caricata, la mappa verrà mostrata all’interno del georeferenziatore. Possiamo ingrandire e rimpicciolire la visuale e muoverci al suo interno con il mouse (e il tasto centrale per spostarsi). La prima cosa da fare è selezionare dei punti di controllo per “agganciare” la mappa al nostro progetto. Per inserire un punto di controllo, dobbiamo innanzitutto assicurarci di aver selezionato il pulsante “Aggiunge un Nuovo Punto” (l’incrocio con un asterisco bianco su sfondo giallo). Per una mappa di questo tipo basterebbero anche solo 6 punti di controllo per georeferirla. Tendenzialmente si comincia da un angolo della mappa e si procede in senso orario (ma questo non è assolutamente una regola).

Georeferire una mappa topografica con coordinate precise

Scegliamo poi un elemento che sia adatto come punto di controllo. In questo caso siamo facilitati dalla griglia di coordinate, che possiamo sfruttare per ottenere i sei punti necessari. Se guardiamo attentamente infatti, le coordinate presentano dei valori. Questi sono valori di northing ed easting, espresse in un sistema di riferimento particolare ma simile all’UTM. Se viene pubblicata una mappa di questo tipo, solitamente viene inserito anche un capitolo dedicato alla descrizione del sistema di riferimento utilizzato sul campo e sulle mappe. In questo caso, la pubblicazione ci informa che la griglia si basa su sulla Israeli Cassini-Soldner (ICS) grid. Il sistema di riferimento di questa griglia è definito come EPSG 28193 (Palestine 1923/Israeli CS Grid), un sistema meno recente e poco accurato ma ancora utilizzato in alcune mappa israeliane.

Questo ci sarà utile da sapere nel momento in cui andremo a completare la georeferenziazione, in quanto dovremo dare a QGIS un sistema di riferimento su cui operare.

Georeferire una mappa topografica con coordinate precise

Cominciamo quindi a cambiare alcune impostazioni (è sempre necessario farlo PRIMA di iniziare a collocare punti di controllo). Apriamo le impostazioni cliccando sull’icona ad ingranaggio in alto (impostazioni di trasformazione).

Georeferire una mappa topografica con coordinate precise

Nel pannello che si apre possiamo modificare il tipo di trasformazione (ovvero come qgis manipolerà la nostra immagine), il sistema di riferimento con cui vogliamo georeferire la nostra immagine, il nome del file di destinazione e il tipo di ricampionamento. Ci sono poi altre impostazioni che al momento non ci interessano ma che potete approfondire nel link alla fine della presentazione.

Georeferire una mappa topografica con coordinate precise

Cominciamo cambiando il tipo di trasformazione. Apriamo il menu a tendina a destra di “Tipo di Trasformazione” e selezioniamo “Polinomiale 2”. Troverete una lista dei modi di trasformazione sempre nei link alla fine, ma per adesso diciamo solo che “Polinomiale 2” è quello suggerito per la georeferenziazione di mappe (anche Lineare sarebbe andato bene in questo caso).

Georeferire una mappa topografica con coordinate precise

Cambiamo poi il sistema di riferimento stavolta sul globo a destra di “SR di destinazione” e nel menu che si apre filtriamo cercando “28193”. Selezioniamo quindi “Palestine 1923 / Israel Cassini Grid” e clicchiamo su OK.

Georeferire una mappa topografica con coordinate precise

Per il momento lasciamo il resto delle impostazioni così e clicchiamo su Ok. Notate che grazie alle due spunte in basso QGIS salverà i nostri punti di controllo e caricherà il risultato una volta eseguito.

Georeferire una mappa topografica con coordinate precise

Possiamo adesso cominciare a inserire i nostri punti di controllo. Clicchiamo ad esempio sull’intersezione delle linee della griglia in alto a sinistra nella mappa , tra 182 e 180. Cercate di essere abbastanza precisi anche ingrandendo l’immagine fino a sgranare un pochino per vedere correttamente l’intersezione.

Georeferire una mappa topografica con coordinate precise

Nel pannello che si apre, potremo decidere di inserire manualmente le coordinate nel campo X e Y, oppure di selezionare un punto sulla mappa. In questo caso, abbiamo già le coordinate a disposizione. Per convenzione, le coordinate riportate sulla mappa sono abbreviate. Se vi ricordate infatti, in un sistema UTM le coordinate Est (Easting) possiedono 6 numeri, e le coordinate nord (Northing) ne possiedono sette. Dato che non abbiamo altre informazioni, diamo per scontato che i restanti 3-4 numeri siano tutti zero. Un’altra convenzione tipica di questo sistema di riferimento è che le coordinate Nord (a destra e a sinistra della nostra mappa) perdono un “1” iniziale quando vengono trascritte. Questo vuol dire che di fronte ad esempio a “180”, dovremo aggiungere un 1 per ottenere il valore corretto delle coordinate.

Georeferire una mappa topografica con coordinate precise

Cominciamo quindi a scrivere la coppia di coordinate, nel campo “X/Est” scriviamo “182000” e nel campo “Y/Nord” scriviamo “1180000” e clicchiamo su Ok.

Georeferire una mappa topografica con coordinate precise

Vedete che il primo punto è stato aggiunto ed appare anche nella lista in basso nel nostro georeferenziatore. Le coordinate che abbiamo appena inserito sono visualizzabili sotto “Dest.X” e “Dest.Y” ovvero le coordinate di destinazione.

Georeferire una mappa topografica con coordinate precise

Aggiungiamo adesso altri 6 punti per completare il quadro (7 punti è il minimo che permette alla trasformazione polinomiale 2 di funzionare), ricordandoci degli zeri da aggiungere alla fine di ogni coordinata X e Y e dell’uno da aggiungere all’inizio di ogni coordinata Y. Prendiamo ad esempio un punto in alto al centro, un punto in alto a destra, un punto in basso a destra, un punto in basso al centro, un punto in basso a sinistra e un punto al centro della mappa. Ricordatevi di seguire un senso orario e di coprire l’estensione della mappa e non solo una parte o un angolo.

Georeferire una mappa topografica con coordinate precise

Una volta completato, vedrete che altre parti della tabella in basso si sono popolate. Quello che ci interessa qui sono i campi “dX”, “dY”, “Residuo (pixel)” nella tabella e il valore “Errore medio” in basso a destra

Questi valori ci indicano eventuali errori nella registrazione dei punti di controllo. I primi due campi (dX e dY) indicano la distanza tra la coordinata che abbiamo assegnato durante la registrazione dei punti e dove quel punto andrà a cadere una volta completato il processo di georeferenziazione. Il campo “Residuo (pixel)” è l’espressione in pixel di questa differenza. L’errore medio è un calcolo matematico rappresentativo di tutte queste differenze. Un valore di errore medio alto, e un valore molto alto nel campo “Residuo (pixel)” indica di solito un errore umano nella registrazione dei punti di controllo, che andranno quindi corretti prima di georeferire l’immagine. Potremo ad esempio spostare un punto di controllo, aggiungerne altri, o rimuoverne alcuni per vedere se la nostra situazione migliora.

Georeferire una mappa topografica con coordinate precise

Sebbene non esista un valore preciso sotto il quale una georeferenziazione è valida o meno (dipende da tanti fattori), tendenzialmente è buona pratica tenere l’errore Residuo (pixel) sotto i 3-4 pixels, ma a volte anche sotto i 10 è accettabile. Un errore alto porta a spostamenti sensibili della mappa oltre che a forti distorsioni.

In questo caso, essendo la griglia sufficientemente precisa e la nostra trasformazione potente, l’errore dovrebbe essere molto basso. Se cambiate trasformazione, anche l’errore cambierà.

Georeferire una mappa topografica con coordinate precise

Possiamo adesso completare la nostra georeferenziazione cliccando sul pulsante con freccia verde (notate come i gcp sono anche visibili sulla mappa).

Georeferire una mappa topografica con coordinate precise

Il processo dovrebbe essere veloce e la nostra mappa verrà caricata nel progetto. Come vedete il risultato è più che accettabile. Possiamo quindi ritenerci soddisfatti e chiudere il georeferenziatore.

Georeferire una mappa topografica non precisa

Georeferire una mappa topografica non precisa

Tuttavia, non sempre sarà così semplice georeferire una mappa. Vediamo adesso un altro esempio con una seconda mappa topografica. Riapriamo la tabella attributi del layer bbox_mappe e scegliamo (cliccando con il numero corrispondente alla riga della tabella) la geometria con nome Nigin e clicchiamo su Zoom alla selezione (il pulsante con la lente d’ingrandimento).

Questo sposterà la visuale della nostra mappa sulla zona coperta dalla geometria selezionata.

Georeferire una mappa topografica non precisa

Una volta fatto, ci troveremo in bassa mesopotamia, nell’area dei siti di Nigin, Lagash e Girsu. Procederemo a georeferire una mappa topografica, per fare ciò, avremo bisogno di poter vedere chiaramente strade, ponti e centri abitati, che però su un’immagine satellitare come quelle offerte da Google, tendono a non essere chiaramente visibili, soprattutto a scale molto grandi di osservazione.

Per risolvere, carichiamo un layer OpenStreetMap andando su Web → QuickMapServices → OSM → OSM Standard

OpenStreetMap è un progetto collaborativo per la creazione di mappature globali a contribuzione aperta a tutti gli utenti. Oltre all’importanza del progetto a livello mondiale, è molto utile in zone dove, ad esempio, Google Maps non è aggiornato a sufficienza.

Georeferire una mappa topografica non precisa

Una volta caricato il nostro layer openstreetmap, procediamo a georeferire l’immagine. Apriamo quindi il georeferenziatore e carichiamo la nostra immagine utilizzando il pulsante “Aggiungi raster” in alto a sinistra.

Georeferire una mappa topografica non precisa

Selezioniamo il file denominato h-38-05, una mappa topografica sovietica redatta intorno agli anni ’80 e scaricabile dal sito https://maps.vlasenko.net/. Una volta selezionata, clicchiamo su Apri.

Georeferire una mappa topografica non precisa

Cambiamo prima di tutto le impostazioni di georeferenziazione, apriamo il pannello delle impostazioni.

Georeferire una mappa topografica non precisa

Qui clicchiamo su “Tipo di trasformazione” e selezioniamo “Helmert”.

Georeferire una mappa topografica non precisa

Cambiamo poi il sistema di riferimento, clicchiamo sul globo e nel filtro scriviamo “32638” e selezioniamo “WGS 84 /UTM Zone 38N EPSG:32638”.

Georeferire una mappa topografica non precisa

Lasciamo tutte le altre impostazioni di default e chiudiamo le impostazioni

Georeferire una mappa topografica non precisa

Assicuriamoci di aver selezionato il pulsante “Aggiunge un Nuovo Punto” (l’incrocio con un asterisco bianco su sfondo giallo) e iniziamo a cercare dei punti di controllo. Sebbene la mappa presenti anche qui una griglia di coordinate, il sistema di riferimento usato (Pulkovo 1942) è di difficile utilizzo e spesso comporta forti distorsioni. Per questo motivo cercheremo sulla mappa degli poi un elemento che sia adatto come punto di controllo. Come detto in precedenza, ponti, incroci e simili sono ottimi punti di controllo in una mappa a questa scala.

Possiamo ad esempio vedere un incrocio ben chiaro in basso a sinistra nell’immagine. Zoomiamo (con la rotella del mouse) in quella direzione.

Orientiamoci e zoomiamo anche nella nostra mappa di base per raggiungere la stessa zona.

Georeferire una mappa topografica non precisa

Orientiamoci sulla mappa di base OpenStreetMap per cercare lo stesso punto, una volta trovato, clicchiamo con il tasto sinistro del mouse al centro dell’incrocio e si aprirà una nuova schermata. Qui, se fossimo a conoscenza delle coordinate del punto, le potremmo inserire manualmente, tuttavia, non essendone a conoscenza, clicchiamo sul pulsante “Dalla Mappa”, ricordandoci di togliere la spunta dalla casella “Nascondi Automaticamente la finestra del georeferenziatore”. Assicuriamoci anche che il sistema di riferimento del punto corrisponda a quello del progetto.

Georeferire una mappa topografica non precisa

Clicchiamo sul centro dell’incrocio anche sulla mappa di base del nostro progetto, e vedremo la finestra di Inserisci Coordinate popolarsi con le coordinate del punto scelto. Se siamo soddisfatti, clicchiamo su ok, altrimenti clicchiamo nuovamente “Dalla Mappa” per correggere eventuali errori.

Georeferire una mappa topografica non precisa

Il nostro punto verrà aggiunto alla lista dei punti di controllo nella Tabella GCP in basso. Selezioniamo adesso un secondo punto. In alto a sinistra nella mappa si trova un centro abitato, dove sicuramente, essendo nei pressi di un fiume, si troverà un ponte. Muoviamoci in quell’area con la mappa del georeferenziatore e facciamo lo stesso nella mappa di base OpenStreetMap.

Georeferire una mappa topografica non precisa

Come pensavamo, un ponte è effettivamente presente sia nella mappa sovietica che nella mappa moderna. Clicchiamo al centro del ponte sulla mappa nel georeferenziatore e selezioniamo sempre “Dalla Mappa” nella finestra di Inserisci Coordinate Mappa. Clicchiamo al centro del ponte anche nella nostra mappa di base e poi clicchiamo su Ok per confermare la nostra scelta.

Georeferire una mappa topografica non precisa

Cerchiamo poi un terzo punto, la parte a destra della mappa potrebbe essere difficile da leggere, ma in alto quasi al centro possiamo vedere chiaramente un altro incrocio, cerchiamo di identificarlo anche su OpenStreetMap seguendo la strada principale verso Est.

Georeferire una mappa topografica non precisa

Individuato l’incrocio (su OpenStreetMap alcune strade potrebbero diventare visibili solo zoomando un po’) ripetiamo il processo per aggiungere un punto dalla mappa e una volta soddisfatti, clicchiamo su Ok.

Georeferire una mappa topografica non precisa

Cerchiamo di localizzare un quarto punto. A destra dell’immagine è possibile individuare un ulteriore centro abitato nei pressi di un canale. Se zoomiamo, troveremo anche qui un ponte utilizzabile come punto di controllo. Individuiamo ancora una volta il luogo anche su OpenStreetMap e zoomiamo per identificare un punto adeguato.

Georeferire una mappa topografica non precisa

Possiamo utilizzare uno dei due ponti presenti nel centro abitato, è indifferente. Ripetiamo ancora una volta il processo per scegliere il punto di controllo dalla mappa e clicchiamo su Ok.

La mappa potrebbe già essere georeferita così, tuttavia, potete cercare altri due o tre punti da utilizzare nella mappa per ottenere un risultato migliore.

Georeferire una mappa topografica non precisa

Una volta cambiate le impostazioni e scelto i punti, possiamo georeferire la nostra immagine cliccando il pulsante verde in alto. In questo caso l’errore è più alto rispetto alla mappa precedente (abbiamo anche un numero minore di punti e una trasformazione meno potente). Se volete modificare i punti, spostarli o rimuoverne qualcuno, lo possiamo fare con i rispettivi pulsanti in alto.

Georeferire una mappa topografica non precisa

Possiamo vedere che la nostra mappa è stata caricata in QGIS e posizionata correttamente secondo i punti che abbiamo scelto fino ad ora. Volendo, se non siamo soddisfatti del risultato, potremo sempre aggiungere nuovi punti e far ripartire la georeferenziazione fino ad ottenere il risultato che vogliamo. Teniamo ovviamente presente che, a causa della scala della mappa e della quantità dei punti di controllo, ci sarà sempre un errore intrinseco nella mappa georeferita, che possiamo cercare di ridurre il più possibile con quanto detto sopra. Chiudiamo pure il georeferenziatore.

Georeferire una mappa topografica non precisa

Possiamo poi eliminare anche i bordi neri dell’immagine dal pannello di stile del layer, scegliendo la scheda della trasparenza e impostando 0 nel campo “Valori nulli aggiuntivi”. Fate attenzione che questo eliminerà anche tutte le scritte in nero presenti nella mappa.

Georeferire una mappa geomorfologica usando punti già esistenti

Georeferire una mappa geomorfologica usando punti già esistenti

Prepariamoci adesso a georeferire la nostra seconda mappa, una mappa geomorfologica. Apriamo nuovamente la tabella attributi del layer bbox_mappe e selezioniamo Sweyhat e nuovamente zoomiamo alla selezione. Ci troviamo nella zona del lago Assad, e la mappa su cui lavoreremo proviene dagli studi della missione americana nell’area di Tell es-Sweyhat, importante sito del III millennio (Wilkinson 2004).

Georeferire una mappa geomorfologica usando punti già esistenti

Facciamo un po’ di spazio per la nostra immagine togliendo la visibilità a tutti i layer che abbiamo caricato in QGIS. Per farlo clicchiamo sull’icona a forma di occhio nel pannello dei layer e dal menu a tendina selezioniamo “Nascondi tutti i layer”.

Georeferire una mappa geomorfologica usando punti già esistenti

Apriamo nuovamente il georeferenziatore (Raster → Georeferenziatore), clicchiamo sul pulsante “Apri Raster” e selezioniamo il file wilkinson_sweyhat.jpg, cliccando poi su Apri.

Georeferire una mappa geomorfologica usando punti già esistenti

Questa volta, invece di inserire manualmente i punti di controllo, li caricheremo da un file. Immaginiamo che qualcuno abbia già georeferito l’immagine, ma che non trovi più il file georeferito. Avere un file con i punti di controllo è molto utile per riprodurre una georeferenziazione in modo rapido.

Clicchiamo sul pulsante a forma di incrocio con una freccia verso l’alto (“Carica Punti GCP”).

Georeferire una mappa geomorfologica usando punti già esistenti

Selezioniamo il file “wilkinson_sweyhat.jpg. points” (nella cartella raster del nostro progetto)

Georeferire una mappa geomorfologica usando punti già esistenti

Una volta caricato, potete vedere che alcuni punti sono stati caricati sulla mappa, possiamo poi osservare la corrispondenza tra questi e il layer dei siti archeologici nel nostro database. Alcuni di questi siti infatti sono presenti nel nostro geopackage e possono aiutarci a georeferire la mappa. Cambiamo intanto le impostazioni del georeferenziatore. Helmert dovrebbe essere già selezionata come trasformazione, quindi cambiamo solamente il sistema di riferimento in “EPSG:32638 - WGS 84 / UTM zone 38N” che dovrebbe essere già visibile dal semplice menu a tendina.

Georeferire una mappa geomorfologica usando punti già esistenti

Come potete vedere, in questo caso l’errore di ogni singolo punto è molto alto (come anche evidenziato dalle linee rosse che partono dai punti nella mappa del georeferenziatore). In questo caso, data la tipologia di mappa potrebbe essere più complicato trovare più punti, quindi al momento accontentiamoci e avviamo la georeferenziazione.

Georeferire una mappa geomorfologica usando punti già esistenti

Come potete vedere dal risultato, la nostra mappa è stata georeferita correttamente, sebbene ci siano delle incongruenza tra i siti nel nostro geopackage e quelli nella mappa, dovute appunto agli errori residui. Potremo riprovare a prendere più punti, oppure a cambiare tipologia di trasformazione (a Polinomiale 1 o 2, ad esempio) per ottenere un risultato, ma non lo faremo oggi. Chiudiamo pure il georeferenziatore.

Georeferire una foto satellitare

Georeferire una foto satellitare

Apriamo nuovamente la tabella attributi del layer bbox_mappe e zoomiamo questa volta nell’ultima area, quella di Tiro. Per quest’area, fortemente e lungamente urbanizzata, non sono disponibili mappe archeologiche o geomorfologiche, ma possiamo ricorrere ad un altro tipo di fonti: le foto satellitari.

Georeferire una foto satellitare

Generalmente le immagini CORONA sarebbero la scelta più comune, tuttavia oggi utilizzeremo delle immagini HEXAGON, che hanno succeduto le CORONA alla fine degli anni ’70 e operato a livello globale durante gli anni ’80. Queste immagini hanno una risoluzione spaziale maggiore delle CORONA e permettono come queste di visualizzare il paesaggio in anni precedenti la grande urbanizzazione del Vicino Oriente.

È stato recentemente pubblicato un articolo molto interessante sulle immagini HEXAGON con dei casi-studio nel Vicino Oriente (Hammer et al. (2022)).

Georeferire una foto satellitare

Apriamo nuovamente il georeferenziatore e carichiamo il raster “hexagon_tyre.tif.”

Georeferire una foto satellitare

Potete vedere che QGIS ha già inserito automaticamente i punti di controllo. Questo perché, se esiste un file .points con lo stesso nome dell’immagine che stiamo georeferendo, qgis li caricherà automaticamente nel momento in cui carica anche l’immagine stessa.

Tuttavia, potete vedere come, nonostante stavolta abbiamo a disposizione più punti delle altre due mappe, l’errore medio è comunque alto. Questo è dovuto alla distorsione dell’immagine (che ricordiamo, è una foto fatta dallo spazio), che necessita di un’adeguata trasformazione per essere corretto al meglio.

Georeferire una foto satellitare

Apriamo quindi le impostazioni di georeferenziazione (il pulsante con l’ingranaggio). In Tipo di trasformazione, apriamo il menu a tendina e selezioniamo Polinomiale 1, lasciando poi tutte le altre impostazioni di default. Una volta fatto, clicchiamo su Ok.

Georeferire una foto satellitare

Notate come l’errore residuo di ogni singolo pixel e l’errore medio siano diminuiti considerevolmente. Generalmente, questo tipo di immagini richiedono molti più punti di controllo, ma per il momento possiamo fermarci qui e avviare la georeferenziazione (che potrebbe richiedere del tempo in più rispetto alle precedenti).

Georeferire una foto satellitare

Una volta caricata la mappa potremo visualizzare il territorio in un’istantanea del 1979, permettendoci di tracciare cambiamenti e modifiche del paesaggio nel tempo se confrontate con la nostra mappa Google di base.

Abbiamo quindi visto come, tramite tre semplici tipologie di mappe, sia possibile visualizzare e importare in QGIS molte tipologie di dati diversi. L’immissione di dati di diverso tipo apre le porte a moltissime possibilità di digitalizzazione e analisi, tutto all’interno di QGIS stesso.

Bigliografia

Davis, M. and Kent, A. J. (2018). Identifying Metadata on Soviet Military Maps: An Illustrated Guide. In M. Altić I. J. Demhardt and S. Vervust (eds.), Dissemination of Cartographic Knowledge, Springer International Publishing, Cham.
Franklin, K. and Hammer, E. (2018). Untangling Palimpsest Landscapes in Conflict Zones: A "Remote Survey” in Spin Boldak, Southeast Afghanistan. Journal of Field Archaeology 43: 58–73.
Hammer, E., FitzPatrick, M. and Ur, J. (2022). Succeeding CORONA: Declassified HEXAGON Intelligence Imagery for Archaeological and Historical Research. Antiquity 1–17.
Kent, A. J. and Davies, J. M. (2013). Hot geospatial intelligence from a Cold War: The Soviet military mapping of towns and cities. Cartography and Geographic Information Science 40: 248–253.
Postnikov, A. V. (2002). Maps for Ordinary Consumers versus Maps for the Military: Double Standards of Map Accuracy in Soviet Cartography, 1917-1991. Cartography and Geographic Information Science 29: 243–260.
Rondelli, B., Stride, S. and García-Granero, J. J. (2013). Soviet Military Maps and Archaeological Survey in the Samarkand Region. Journal of Cultural Heritage 14: 270–276.
Wilkinson, T. J. (2004). On the Margins of the Euphrates: Settlement and Land Use at Tell Es-Sweyhat and in the Upper Lake Assad Area, Syria, Oriental Institute Publication 124, Chicago.
Davis, M. and Kent, A. J. (2018). Identifying Metadata on Soviet Military Maps: An Illustrated Guide. In M. Altić I. J. Demhardt and S. Vervust (eds.), Dissemination of Cartographic Knowledge, Springer International Publishing, Cham.
Franklin, K. and Hammer, E. (2018). Untangling Palimpsest Landscapes in Conflict Zones: A "Remote Survey” in Spin Boldak, Southeast Afghanistan. Journal of Field Archaeology 43: 58–73.
Hammer, E., FitzPatrick, M. and Ur, J. (2022). Succeeding CORONA: Declassified HEXAGON Intelligence Imagery for Archaeological and Historical Research. Antiquity 1–17.
Kent, A. J. and Davies, J. M. (2013). Hot geospatial intelligence from a Cold War: The Soviet military mapping of towns and cities. Cartography and Geographic Information Science 40: 248–253.
Postnikov, A. V. (2002). Maps for Ordinary Consumers versus Maps for the Military: Double Standards of Map Accuracy in Soviet Cartography, 1917-1991. Cartography and Geographic Information Science 29: 243–260.
Rondelli, B., Stride, S. and García-Granero, J. J. (2013). Soviet Military Maps and Archaeological Survey in the Samarkand Region. Journal of Cultural Heritage 14: 270–276.
Wilkinson, T. J. (2004). On the Margins of the Euphrates: Settlement and Land Use at Tell Es-Sweyhat and in the Upper Lake Assad Area, Syria, Oriental Institute Publication 124, Chicago.